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Penetration of X- and Gamma Rays to Extremely 

Great Depths' 

U. Fano ' 

Earlier work on the asymptotic trend of the X-ray intensity at great distances from a 
source is reviewed and completed in various aspects. The asymptotic law is shown to be the 
same as in the '^straight-ahead'' approximation (which disregards deflections) whether the 
primary energy is higher or lower than the energy of minimum absorption, provided a constant 
is replaced by the eigenvalue of a suitable Wick equation. The penetration in directions 
oblique to the source direction hardly ever attains its asymptotic trend when the source 
energy is lower than the energy of minimum absorption. This situation raises a difficult 
problem regarding the penetration law in the range of great depths where the asymptotic 
trend is being approached very slowly. 

1. Introduction 

The deep penetration of X-rays, or gamma rays, in an infinite, homogeneous mecHum 
are discussed. Much fundamental understanding of this phenomenon derives from the work 
of Wick [1] ^ on the analogous phenomenon of neutron penetration. Vice versa, the present 
work might contribute additional clarification to neutron prol)h'ms as well as somc^ background 
for still unsolved problems of charged-particle penetration. 

This paper constitutes a final report of developments since 1948. It includes a review^ of 
preliminary reports [2, 3, 4] and new material tliat wasrequircnl to complete^ the initial ])rogram. 
Some concepts and techniques developed in the course of tlu* study but not utilized in the even- 
tual solution are nevertheless reported briefly. 

The assumption of an infinite, homogeneous inedium disrt^gards tlu* eff(»ct of boundaries 
and inhomogeneities, wd)ich constitutes a separate, still larg(*ly unsolved, problem. 

The very deep penetration of X-rays depends primarily on the course of multiple C\Mn])ton 
scattering under conditions where photoelectric absorption and pair production are compara- 
tively unimportant. At the low-energy end of the spec^trum, for example, below 50 kev, wdiere 
the energy shift of the scattered photons can be disregarded or treated as a small correction, 
the diffusion of photons has been studied by Chandrasekhar [5], but w ithout specific reference 
to very deep penetrations. (Notice that photons below 50 kev disappear rapidly by photo- 
electric effect, in most materials.)^ At the high-energy end, where large amounts of X-rays are 
rc^generated by electrons, a complete study of X-ray penetration would require a treatment of 
the whole cascade shower process. The X-ray regeneration remains moderate for secondary 
electrons up to 10 Mev in lead, and up to 100 Mev in light materials, and may be treated, if 
necessary, in this energy range as a secondaiy source. The present paper disr(*gai-ds the re- 
generation of X-rays by electrons and treats pair production as a mechanism of outright absorp- 
tion. Nevertheless, its results apply to the penetration of the tail end of showers that is con- 
trolled by photons below 10 Mev in lead and below 100 Mev in light elements. 

We shall deal primarily with photons from 50 kev to 50 Mev which experience concurrent 
processes of energy degradation and multiple scattering terminating in outright absorption. 
Repeated Compton scattering of the ^^primary" radiation emitted by a source gives rise to 
'^secondary" X-rays of lower energy traveling in all directions. In a medium of low atomic 
weight, a photon may be scattered 5 or 10 times, on the average, before eventual absorption 
by photoelectric effect. The main complication in our study arises from the generation of 
this complex radiation. 



' Work supported by the Office of Naval Research and by the Atomic Energy Commission Reactor Division. 

2 Appendix C is by L. V. Spencer. 

3 Figures in brackets indicate the literature references at the end of this paper. 

^ This consideration applies also to the further development of Chandrasekhar's method that was carried out by O'Rourkc [21j in connection 
with the X-ray penetration problem. 
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This complication and the compHcated dependence of the cross sections on the photon 
energy preclude a fully analytical treatment of the problem. Therefore, the analytical treat- 
ment that is developed in this paper does not lend itself readily to direct numerical applications. 
However, its results provide one of the footings for a practical method of numerical calcula- 
tion [6] whose range of application extends to very large distances from the source. The 
analytical laws of intensity variation versus depth also serve to extrapolate the results of num- 
erical calculations to still greater distances from the source. 

Throughout the present investigation, much reliance has been placed on gaining an 
initial qualitative understanding of the factors that control the process of degradation, scat- 
tering, and penetration of X-rays. Such an understanding enables one to adapt the mathe- 
matical procedure to the conditions affecting each specific problem. In this respect our work 
departs from other work in the same field that has been characterized by the initial adoption 
of less flexible procedures. For example, the X-ray distribution has been resolved into '^orders 
of scatterings^ (that is, into once-, twice-, thrice-, . . ., scattered photons) [7, 8, 9]. The ex- 
pansion converges poorly where the photons experience protracted scattering and its applica- 
tion becomes cumbersome. However, at least in some instances, adequate calculations of 
penetration have been made [8]. Semiempirical improvements to the procedure have also 
been introduced [9]. Alternatively one may schematize the scattering and degradation 
process in such a way that the corresponding basic equations are sufficiently simplified to 
allow analytic solution [7, 9, 10, 11]. These simplifications, however, seem generally to in- 
volve considerable limitations in the applications of the theory and may delete some char- 
acteristic features of the phenomenon. Finally, calculations by the ^'Montecarlo'' sampling 
methods are possible [12] but do not yet appear convincingly successful, at least in their appli- 
cation to very deep penetrations, which requires the use of '^biased sampling' \ 

Our treatment will be limited to X-ray distributions having a plane symmetry, that is, 
which vary in one direction of space only and are generated by sources having the same sym- 
metry. The study of spherically symmetrical distributions is essentially equivalent to that of 
plane-symmetrical ones. The relationship between these types of distributions is well known 
[1] and will be indicated in Appendix A. We shall also assume that the source is concentrated 
on a single plane, which implies no additional restriction owing to the linearity of the problem. 
Notice that plane symmetry may always be attained by considering, instead of the X-ray 
flux at each point and in each direction, only the integral or the average flux over all points 
of a plane and over all directions forming the same angle with the plane. The reason is that 
this plane integral flux generated by a localized source is equal to the localized flux generated 
by a plane distribution of sources. 

The most general type of source may be regarded as an aggregate of point monodirectional 
sources. The study of X-ray penetration under conditions of plane symmetry may serve as a 
first step for the treatment of point monodirectional sources. This approach has been ex- 
ploited in the study of X-ray penetration to small or moderate distances from a source [13]. 

2. Survey of the Problem 

This section contains a qualitative analysis of the deep penetration of X-rays, with the 
purpose of pointing out the specific problems that require a detailed mathematical treatment. 

2.1. Transient Processes and Equilibrium States 

The layers of material near an X-ray source are traversed almost exclusively by primary 
radiation because it takes a certain thickness of material to generate a substantial amount of 
secondary radiation. Therefore, one expects the intensity of secondary radiation to build 
up rapidly from one layer to the next near the source in a sort of transient process. Vice versa, 
at greater depth within a material one may hope to find some sort of steady state in which the 
secondary radiation is present in substantial amounts. The nature of this steady state is the 
main object of this investigation. 
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A typical steady state in wliicli a softcM- I'adiation (•ontinually arises from a liar-(l(M- one is 
commonly called a state of radiative '^ equilibrium." The X-ray penetration problem does not 
lead to a typical state of equilibrium, but it is helpful to review first the properties of such a 
state. Equilibrium involves, in general, three different features, namely: (1) The ratio of the 
intensity of the secondary to the intensity of the primary radiation approaches asymploticalh' 
a maximum value, as the depth of penetration increases. (2) The quality of the scH'ondary 
radiation also becomc^s independent of the depth of penetration, because the intensities of its 
various spectral components ])ear a constant asymptotic ratio to the primary intcvnsity. (3) 
Therefore, the intensity of the whole radiation is controlled b^^ the progressive attenuation of 
the primary radiation, which follows an exponential law in the case of X-rays. Tf these 
circumstances applied, they would greatly simplify the stud^" of penetration. The maxinunn 
values of intensity ratios mentioned in (1) and (2) are determined by the ratios b(^tw(HMi the 
attenuation coefficients ^ jjl of the X-rays of various energies. The depth of penetration at 
which the intensity ratio between any two X-ray components attains a value within, say 10 
percent of the maximum, turns out to be inversely related to the difference between the attenu- 
ation eoefRcients of those components. 

2.2. State of Limited Equilibrium 

The approach to a steady state under the conditions of hard X-ray penetration presents 
the following characteristic departure from tlie more familiar equilibrium conditions. In the 
Compton effect, some of the X-rays have very nearly the same energy, the same direction, and 
the same attenuation coefficient after as before the scattering. In fact, the attenuation coeffi- 
cient of some among the secondary X-ray components differs only injinitesimally from the 
attenuation coefficient of the most penetrating X-rays that are present. This circumstance 
does not mereh^ slow down the approach to equilibrium, it actually sup])]'(^sses tlu^ features (1) 
and (3) in that the intensity ratio of secondar^^ X-ra^' components to a monochromatic primary 
may grow beyond any limit . 

At the same time one may expect that, within a finite^ thickness of shield, some sort of 
equilibrium should be attained among those components of the secondary X-ray spectrum, 
whose attenuation coefficient differs from that of the most penetrating component by a suffi- 
ciently large amount. In fact, the relative intensity of most secondary components and, 
therefore, the average ''quality'^ of most of the secondary radiation turn out to approach a 
limiting value within the thickness of a finite shield, as indicated in feature (2). On the other 
hand, features (1) and (3) do not obtain. The intensity ratio of this secondary radiation to 
the primary X-rays does not approach a maximum value, but it keeps increasing; the attenua- 
tion of the primary radiation alone does not control the attenuation of the secondary radiation. 

In a hypothetical state of true equilibrium, the flow of radiation at great depths within a 
shield is described as the product of an exponential function of the depth and of a fuiu'tion of 
the energy and direction of each secondary radiation component. In the sort of steady state 
that is actually achieved, the flow of most of the secondar}^ radiation at great depths is described 
as a product of a nonexponential function of the depth and of a function of the energy and 
direction of each component. Calling x the depth of penetration, measured from the source, 
E the energy of a photon, and 6 the angle that its direction of travel forms with the x axis, 
Eq the energy of the primary photons, the flux of photons of various energies and directions 
^{x,E,d) takes then the form 

Hx,E,d)^f(x)g(E,d) (1) 

when x is sufficiently large and E sufficiently smaller than Eq. 

At large depths most of the energy is carried by the softer secondary components. One 
can, therefore, use the approximate form of the flux to evaluate the quantities of practical 
interest. These quantities vary, as a function of depth, in proportion to f(x), 

5 The term "attenuation coefficient" is used instead of "absorption coefficient" to make it clearly understood that n includes the effect of 
removal of X-rays from a narrow beam by Compton effect in addition to their true absorption by photoelectric effect or pair production. 

97 



2.3. Intensity Variation at Great Depth 

The over-all attenuation of the radiation does not follow an exponential law because the 
most penetrating secondary components never approach an equilibrium. The determination 
of the over-all course of attenuation, which is described by the function/(ic), requires, therefore, 
a study of the ever developing transient process of formation and destruction of the hardest 
secondaries. This study constitutes the main topic of the present paper. 

By centering one^s attention on the formation and absorption of the most penetrating 
secondaries, one attains a substantial simplification of the problem. The secondaries that 
need be considered cover only a narrow spectral range. Therefore, the variations of the 
scattering and absorption cross sections within this spectral range can be treated as small 
quantities. Furthermore, the greatest contribution to the deep penetration arises from 
secondaries directed in a ver^^ narrow beam perpendicular to the source plane. ^ 

However, the small differences of attenuation coefficient and of direction among the most 
penetrating components cannot be disregarded, even in the first approximation. On the 
contrary, it is just these diferences that control the relative intensity of the various components 
and thereby the over-all course oj attenuation at great depths. The following argument indicates 
that small changes of attenuation coefficient and of direction have comparable importance. 
Radiation of wavelength X, which travels in a direction forming an angle § with the axis of 
penetration, r, travels a distance A/cos ^ as its distance from the source plane increases from 
X to x+A. Therefore, its intensity decays along x with an effective attenuation coefficient 
/x(X)/cos ^. Compton scattering of a photon, with a wavelength shift of 5X Compton units, 
results in a change of both jli(X) and cos ??. If a photon starts with a wavelength Xq in a direction 
with ^0=^0, cos ^0=1, after one scattering its direction changes, according to the Compton 
law, to t? = arccos [1— 5X] and its effective absorption coefficient from /x(Xo)/cost?o=M(Xo) to 
/i(Xo+6X)/[l — 5X]--M(Xo) + (rfAi/^X)o5X+/x(Xo)5X. The corrective term (c//x/rfX)o5X arises from the 
change of /x, the term jLt(Xo)6X from the change of cos ^. The ratio of the two corrective terms 
is d log m/^X, with X in Compton units, which is a number of the order of 1 for most materials 
and wavelengths of interest. 

The importance of the small deflections was not appreciated in the early stages of the 
investigation, which disregarded altogether the changes of direction experienced by the most 
penetrating radiation components, that is, relied on the ''straight-ahead approximation.'^ 
The results obtained in this manner, which are reported in [2, 3] and reviewed in section 5, 
proved useful, nevertheless, because the effect of small deflections does not change the analytical 
form of the intensity variation but only the value of certain numerical constants. This simple 
mode of action of the deflections emerges as a result of the analytical treatment of section 7 
and may be visualized as follows. Photons that travel in directions oblique to the direction 
of maximum penetration are selected-against in the course of deep penetration, much as if 
they had, in effect, a larger attenuation coefficient /i(X) (see above). The selection operates 
simultaneously in favor of photons with near-minimum /x and near-minimum, that is, near-zero, 
obliquit^^ §. As a result, photons of each wavelength are concentrated in a cone of directions 
whose aperture is the narrower the nearer to the minimum is ju. This is to say that the effect 
of small obliquity merely parallels and amplifies the effect of small variations of the attenuation 
coefficient. Indeed, the two effects result in a fixed ratio to one another (see sections 6 and 
7), so that the study of the effect of variations of ii alone, in the straight-ahead approximation, 
gives the correct analytical form for the trend of the total intensity at great depths. 

Two different analytical forms for this trend result under different circumstances, depending 
on the type of variation of the attenuation coefficient among the most penetrating components. 
As long as the energy of the primary X-rays is not too large, none of the secondary X-rays 
is more penetrating than the primary X-rays. The hardest secondaries are then simply those 
whose energy is just a little lower than the energy of the primaries. The progressive accumula- 

• When the whole primary radiation travels in directions oblique to the source plane, there occurs a more complicated situation that will be 
discussed in section 2.6 and in section 8. 
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tion of secondary components under these conditions has been studied in [2] under the 
straight-ahead approximation and in |[4], taking into account the deflection effects. The 
intensity distribution law f{x) is shown there, and in section 9, to be of the form 

f{x) = x^ exp ( ~ Mo^) . (2) 

Here mo is the attenuation coefficient of the primary X-rays and K is a niiinher, usually of the 
order of \ ; K depends on the cross section for Compton scattering and on the first diM'ivative 
of the atti'nuation coefficient with respect to energy, evaluated at the en(»rg3' of the primary 
X-rays. (For numerical values of K and the parameters on which it depends, see appendix 
C, tables 1 and 2.) The exponential factor represents the usual eff'ect of attenuation of the 
primaries, whereas the factor x^ represents the effect of accumulation of hard secondaries and 
is called the ^^buildup factor''. Numerical values of K are given in appendix C. 

At very high energies, absorption by pair production increases so rapidly that higher 
energy X-rays are actually less penetrating than lower energy ones. Therefore, in the case of 
high-energy primaries, some of the secondary X-rays are more penetrating than the primaries. 
The attenuation of the total radiation depends upon the formation and decay of the most 
penetrating secondary components, whos(^ photon energies nuiy be nuich lower than that of 
the primaries (in lead the hardest photons have about 3-Mev energy). This phenomenon has 
been treated in [3] in the straight -ahead a])proximation, })ut the eff'ect of small deflections 
has not been taken into account in previous papers. Both trtuitments are given in sections 5, 
7, and 9 and show the intensity distribution law f{x) to be of the form 

f{x) = x-'^' exp[//(M.4^]oxp(-M.a:). (3) 

Here Hm indicates the mininnun value of the attenuation coefficient, that is, the attenuation 
coefficient of the most penetrating X-rays in the material under consideration. The factor 
exp(— Mm^) represents the attenuation of these X-rays. The factor x~^^^ cx'p[H(fi,nX'] is the 
buildup factor, corresponding to x^ in (2), and His a constant that depends on the cross section 
for Compton scattering and on the values of the attenuation coefficient and of its second 
derivative, evaluated at the energy of the most penetrating X-rays. (For values of M/« and //, 
see appendix C, table 3.) 

2.4. Numerical Methods Combined With Analytical Results 

The analytical results, (2) and (3), represent limiting forms of the intensity distribution, 
valid at extremely great depths, where the distribution depends on the generation and attenua- 
tion of an extremely narrow range of secondary spectral components. Those* results must be 
supplemented in two directions. In the first place, one must find not only the relative X-ray 
intensity at different positions but the absolute intensity, particularly of lower energy compo- 
nents traveling in various directions, that is, the function.^ (£',?>) of eq (1). In the second place, 
one must be able to carry out calculations without excessively stringent limitations on the 
important range of secondary spectral components. 

Various partially successful attempts were made to develop methods of successive approxi- 
mation to take into account successively higher derivatives of the cross sections in the important 
spectral range. However, much greater success was eventually achieved by the semiasymptotic 
numerical method of Spencer [6], which fulfills both of the requirements indicated above. 
Additional numerical results for the X-ray distribution at moderate depths are provided by 
polynomial method of calculation [13]. 

Reliance on numerical methods to the extent of solving the transport equation numerically 
was dictated by the complicated dependence of the cross sections upon the X-ray energy. At 
the same time, since the transport equation involves three independent variables (namely 
distances from the source, photon energy, and direction of propagation), the reduction of the 
numerical burden to manageable proportion requires much guidance from qualitative analysis, 
and as much help from analytical development as^conveniently possible. Spencer has stressed 



the rapid gain in computational efficiency that can be derived from general information on the 
desired solution [6, 14]. Economy also arises from suitable substitutions of the variables of 
the problem (see, for example [13]). 

The space variable, depth of penetration from the source, can be separated from the other 
variables by a Fourier-Laplace transform. This transformation has the additional advantage 
of replacing the depth of penetration, which varies over a wide range, with a transform variable 
whose important range is narrow. As will be seen in section 4, the analytical treatment that 
determines the intcDsity distribution laws (2) and (3) deals specifically with the behavior 
of the distribution function at a singuk- point of the transform variable. The numerical 
method holds only for values of the transform variable at a finite distance from the singular 
point. The analytical treatment serves as an essential complement to determine the critical 
behavior at infinitesimal distances from the singularity. Furthermore, the knowledge of this 
analytical behavior and of the resulting behavior of the corresponding inverse transform, 
(the function /(x) of eq (1), (2), and (3)), serves as a guide for the inverse transform procedure 
that is also required in the numerical work. 

The two remaining variables, photon energy and direction of propagation, can be separated 
only within the limits of application of the analytical treatment, that is, for extremely deep 
penetration and narrow range of photon energies (section 7). Otherwise, the interlinkage of 
energy and direction offers the most serious difficulty to the solution of the problem. The 
difficulty has been overcome in the numerical work by learning how to take into account this 
interlinkage without excessive complication [6]. This is done by describing the directional 
distribution adequately by means of a few parameters; namely, moments when the distribution 
is peaked, and coefficients of a Legendre expansion when it is flat. The interlinkage of energy 
and direction is reviewed briefly in section 6. Numerical methods have to complement the 
analytical treatment even when the separation of variables succeeds in order to work out the 
directional distribution (section 7 and appendix B). 

2.5. Verification of the Qualitative Analysis 

The introductor}" qualitative picture of the X-ray penetration, which has led to the results 
indicated by eq (1), (2) and (3), has been verified to a considerable extent by experimentation 
and by independent numerical calculations. 

Two experiments dealt with a point source of Co^%-rays surrounded by a mass of water. 
In the first experiment [15] ionization-chamber and Geiger-counter measurements of total 
intensity were made at various distances from the source. The results of these measurements 

should be represented, in terms of the approximate eq (1), by /(a:) 27r sin Md^ g(E,i})r(E)dE, 

Jo Jo 

where r(£') indicates the response of the measuring instrument. In the absence of information 
on g{Ejd), the plot of the measurement versus the depth of penetration x should follow the trend of 
the f unction /(x) . The plot does, in fact, take up the trend of (2) at great distance from the 
source (1.5 to 2.5 m of water), with a value of K equal to that predicted by the detailed theory of 
[4] and of section 7. 

In the second experiment [16], the spectrum of the secondary electrons in the water was 
measured at various distances from the source. The plot of the number of electrons of various 

energies e, should be represented according to eq (1) hj J{x) 27rsin??c??> g{E,i})R(E,€)dE, 

Jo Jo 

where R(E,€) indicates the probability of production of an electron of energy e by a photon of 
energy E. Qualitatively, the shape of the plot becomes independent of the penetration depth 
X at great depths as predicted by (1). 

Equation (1) gives no information about the shape or the absolute magnitude of g{E,§). 
However, rather detailed calculations of the complete distribution function ^{x,E,d), for all 
depths of penetration covered by the experiments can be made by the independent method of 
•polynomial expansion [13]. These calculations have successfully predicted the quantitative 
results of the two experiments. 
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FicuRE 1. (a) Differential spectra I{E,iior) of the 
X-ray intensity at various distances r from a Co-60 
point isotropic source in water. 

The higher energy portions of the spectra have been omitted for clarity. 
The parallelism of the curves on the left side of this semilogarithmic plot 
indicates the trend towards equilibrium. 

(b) Directional distribution I{E,^rfr) of the X-ray 
intensity for photons of several energies E at several 
distances r from a point isotropic source of Co-60 in 
water. 
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FiGURE 2. Differential spectra Yo(fioX,E) 
at various distances x from a plane mono- 
directional source of 5.11 Mev photons in 
Fe. 

Notice the trend to parallelism on the left side as in 
figure 1, a. 

3 101 




2 3 

ENERGY, Mev 



In addition, the results of these calculations for moderately great depths provide a semi- 
quantitative check of the predictions of the theory of penetration to very great depths. These 
predictions find good verification. In particular, the trend toward an equilibrium distribution 
for the lower-energy secondary X-rays is clearly displa^-ed both with respect to the spectral [13] 
and to the directional [17] distributions (see fig. 1). The emergence of an equilibrium state at 
lower energies, while the higher energy components persist in a state of nonequilibrium is 
most clearly displa3^ed in the results of semiasymptotic calculations [6] (see fig. 2). 

2.6. Influence of the Source Direction on the Penetration Law 

Previous theoretical work on the penetration of radiation from plane sources has dealt 
onl}^^ with isotropic sources and with sources concentrated in a direction perpendicular to the 
source plane. In general, one may want to deal with source radiation aimed in directions 
that form a fixed skew angle with the source plane. However, this new problem raises questions 
that are largely unsolved and that will be indicated briefly here and discussed somewhat further 
in section 8. 

At the outset one must distinguish between the conditions where the energy of the primary 
X-rays lies below and above the energy of minimum attenuation (corresponding to eq (2) and 
(3), respectively). In the first event the penetration law is controlled by secondary X-rays of 
energies just below the source energy. The overwhelming majority of these secondary X-rays 
travel in directions very nearly as oblique as those of the primary X-rays. Therefore, the 
obliquity of the source has a very major influence on the intensity distribution at very great 
depths. 

In the second case, when the source energy exceeds the energy of the most penetrating 
secondaries by a substantial amount, the obliquity of these secondaries is less immediately 
related to the obliquity of the primaries. It will be shown in section 8 that the penetration is 
controlled by the accumulation of secondaries that have been degraded by steps to the energy 
of minimum attenuation and simultaneously deflected to directions nearly perpendicular to 
the source plane. The absolute intensity of these secondaries depends, of course, on the 
obliquity of the source but their rate of accumulation and, hence, the law of deep penetration 
are independent thereof. 

Therefore, the main difficulties arising from source obliquity occur only when the energy 
of the primaries lies below the energy of minimum a tteu nation. If ?^o indicates the angle 
formed by the source direction with the x direction, perpendicular to the source plane, the 
*^ effective attenuation coeflftcient" of the primary X-rays of wave length \ is m(^o)/cos ?^o (see 
section 2.3). After one Compton scattering with a wavelength shift d\, the obliquity ?> is 
given by cos ??=cos ??o[l — 5X] + sin ?^o[25X— 5X^]^^^cos (p, where ^ may have any value. Minimum 
effective attenuation after one scattering is obtained when <p=0, but this minimum value 
Ai(Xo+5X)/<cos??>maxniay be larger or smaller than the initial value jLt(Xo)/cos i^o, depending on 
the values of djji/dX and of ^^o- The most favorable reduction of the effective attenuation co- 
eflftcient is brought about by a succession of scatterings with infinitesimal 5X^s, all directed so 
as to reduce ^. An infinite sequence of such scatterings would eventually bring about a deflec- 
tion all the way to ^^0 without any finite increase of /jl, but the probability of such an event 
is, of course, infinitesimal of a high order. Nevertheless, an unlikely process of this kind must 
determine the penetration at extremely great depths. That is the penetration will be eventually 
controlled by the scanty accumulation of secondary X-rays with energies very near the pri- 
maries and direction nearly perpendicular to the source plane, that is, with effective attenuation 
coefficients near jjlq. The course of accumulation of secondary X-rays that are effectively more 
penetrating than the primaries because of reduced obliquity has not yet been studied. Knowl- 
edge of the limiting trend at extremely great depths may not be very important because it may 
not even become recognizable except under unrealistic conditions. Nevertheless, lack of 
understanding of this process may undermine the effectiveness of the semiasymptotic and 
polynomial methods of calculation, whose application is now being planned. 
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If the source emits X-rays of various obliquities, some of which travel initially iu the 
direction perpendicular to the source plane, it will be just the X-rays in this direction and 
their secondaries that achieve great penetration. Deeper and deeper penetrations arise from 
narrower and narrower })undles of primary X-rays. Therefore, a source with a broad distribu- 
tion of directions contributes less and less toward deep penetration, as compared to a source 
aimed exactly in the direction perpendicular to the source (?^o=0)- This argument will be 
confirnu^d by the analytical trcvatment of section 8. In particular an isotro])ic source, or any 
source with practically uniform distribution in directions nearly perpenducilar to the source 
plane, yields an intensity distribution law in accordance with eq (2) but with a value of K 
lower by 1 than the value for a source concentrated at ??o=0 [6]. 

3. Transport Equation 

The following symbols w411 be used: 

r=(x,y,z) = coordinates of a point of a medium. 
X= wavelength of a photon in Compton units (h/rnc). 
u=(Ux,Uy,Uz) = {i^,(p)=umt vector indicating a photon direction. 
/x(X) = total (narrow-beam) attenuation coefhcient of the medium for |)h()lons of wave- 
length X. 
Mt/2 = Thomson scattering coefficient (probability per unit path) of the medium^ 
low energy limit of the integral Klein-Nishina coefficient. 
/:(X^X) = (:V8)MT/,[X7X+X'V^H(X7X)^XX-X'-2)(X-X0] = Klein-Nishina differential 

scattering coefficient of the medium for Compton scatt(M"ing with a wave-length 
change from X' to X, per unit X, for X'<X<X' + 2, 
y(r,M,X)= number of photons per unit volume, per unit solid angle and per unit X. 
6(Q = Dirac^s delta function. 
Xo= wavelength of the source photons (if monoclu-omatic). 
/xo=m(Xo). 

X^^wavelength of the photons whose attenuation coc^llicient has the lowest A^alue 
in the range of integration from Xo to X. 

Ms=m(Xs). 

X,^= wavelength of the photons whose attenuation coefficient is the absolute mini- 
mum in the medium under consideration. 

The degradation, penetration, and (lifl'iision of X-rays is governed by a transport equation. 
This equation represents the rate of change of the density of photons with a specified direction 
and wavelength, from one point to the next in the direction of propagation, as the sum of three 
terms namely: (a) the attenuation of the photon density as a result of absorption and scattering, 
(b) the addition of photons that take up the specified direction and wavelength as a result of 
Compton scattering, (c) the addition of photons wntli the specified direction and wavelength 
from the source. The equation is: 

wgradY{r,u,\) = ~iJL{X)Y(r,u,\) 

+Jx''_9^^'^(^''^)j4,^^' (27r)-^5(l-M-a'-X+X0 Y(r,u\y)+S{r,u,\), (4) 

where S(r,u, X) represents the source, that is, the density of photons produced at r with the 
direction u and the wavelength X. The 5-function represents the Compton law which requires 

X to equal X' + (l — mm')- 

Following the plan indicated in the introduction, we assume that the source distribution 
has a plane symmetry, that is, that it does not depend on y and z, nor on the component (p of u, 
(but only on i^ = arccos Ux). As a result, Y will be similarly independent of y,z and (p. We also 
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assume that the source is concentrated on the plane x = and at a wavelength X=Xo, but this 
constitutes no real restriction owing to the linearity of the equation. Accordingly, we have 
S(r,u,\) = 8(x)d(\—\Q)f(Ux) and the equation reduces to 

U:c ^ Y{x,U:cy X) / dx = — jLt ( X) Y(:x,u^, X) 

Notice that the kernel k(y,\) of this equation may be multiplied by any ratio ^(X')/^(X), 
provided Y is suitably renormalized. 

4. Fourier-Laplace Transform 

The transform method involves a separation of variables. It represents the X-ray distri- 
bution as a superposition of components, each of which has a certain angular and spectral dis- 
tribution uniform over all the space. The intensity of each component varies exponentially or 
sinusoidally from point to point according to exp(—px), where p may be complex.^ 

The transform and its inverse are: 

y{p,u^,\) =J_"^ exp ijpx) Y{x,u^,\)dx, (6) 

Y{x,u^,\) = {2Tri)-^ \ .'" exp (—px)y(p,u,,\)dp. (7) 

J — zoo 

K The transform method of analysis consists of two steps: (a) A study of the directional and 
spectral distribution y{p,Ux,\) for various values of ^, (b) an evaluation of the inverse transform 
integral (7). 

The evaluation of the integral may proceed by the path of steepest descent. Whether this 
method is eventually followed or not, an analysis of the ^ topography" of the transform y(p,Ux,X), 
that is, of the distribution of its absolute values over the complex plane p proves very useful. 
This analysis indicates what region of the plane yields the largest contribution to the integral 
when the steepest descent method is followed. The determination of y{p,Ux,\) in this region, 
which is often very limited, proves anyhow to be very efficient for evaluating the inverse 
transform. 

The equation that governs the directional and spectral distribution y{p,Uxj\), for each 
value of p, derives from (5) by means of the transformation (6), that is, by multiplication by 
exp (pa::) followed by integration. This equation is: 

[yi{\)-pux]y{p,u,,\)= rdyk{y,\) f du\2T)-%\-wu~\-\-V)y{p,u^,V)+b{\-\)j{Ux). (8) 

Jo J47r 

The source term in this equation does not depend on p because the source term in (5) is con- 
centrated on the plane x=0. 

Because the wavelength of a photon increases as a result of Compton scattering, the dis- 
tribution of photons of each wavelength X depends on the distribution of photons of shorter 
wavelengths X' but not on that of longer wavelength photons. If one knows the photon dis- 
tribution (in spectrum and direction) at all wavelengths X'<^X, one can evaluate the integral 
in (8) and thereby find the distribution y(p,Ux, X) on the left side of the equation. Thus eq (8) 
may be solved in principle stepwise, proceeding from the source wavelength Xq to longer and 
longer wavelengths. This procedure, which reflects the course of the physical process of energy 
degradation, points to an important effect. If the intensity of X-rays of a certain wavelength 
builds up to a high value, Compton scattering of these X-rays builds up, in turn, the intensity 
of all longer wavelength radiation. 



7 Notice that p has here a sign opposite to that which is most common in the literature. The convention adopted here proves convenient in 
problems of penetration and straggling. 
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Figure 3. Elements of a relief map diagram of the trans- 
form distribution y on the plane of the complex variable p. 
, Lines over which y(p) is singular; — 



trend of the steepest 
descent line of y(p) exp (—px) for positive p;)i, approximate position 
of the saddle-point. 




Figure 4. Diagram of a plot of \y(p) \ along the real axis 
of figure 3. 



Equation (8) shows that y{p,Ux, X) becomes infinitely large when ii{\) — pUx vanishes. Fur- 
thermore, y grows rapidly for those combinations of p,Ux and X for which yL{\) — pUx is particu- 
larly small. These combinations of variables are critical in various respects. We just men- 
tioned that once y{p,Ux, X) has become very large for certain values of u^ and X, it remains 
correspondingly large for all longer wavelengths. We also know that the distribution of large 
values of y(p,Ux, X) determines the path of steepest descent of the integral (7). Finally, the 
process of numerical integration becomes increasingly difficult as the rate of rise of y{p,Ux, X) 
increases, so that one may want to utilize, at least locally, analytical methods of solution. 

The coefficient ^{\) — pUx of (8) vanishes only for positive real values of p larger than /x(X) 
or for negative real values smaller than — /x(X). Therefore, eq (8) beconu's singular only for 
values of the complex variable p which are confined to the sections of the real axis marked 
heavily in figure 3. The lines terminate at p^dz/x^j where Hs indicates the smallest value of 
/i(X') in the range of wavelengths from Xq to X. There appears to be no other point in the plane 
p where y{p^Ux, X) diverges. 

Figure 4 shows a schematic plot of \y\ against real values of p. The plot rises to infinity 
at p= ijUs, where the equation has become singular for Ux=\ and for the w^avelength X5<X at 
which /x takes its smallest value ^x^. 

The integrand of (7) equals y{p,Ux, X) times exp{—px). The factor exp ( — pj::) , regarded as 
a function of p, slopes down to the right for positive values of x and to the left foi- negative 
values. The slope increases in proportion to the distance from the source, \x\. Therefore 
the integrand has a minimum, for real p, where d log y(p,Ux\)/c>p=x, that is, at a point that 
approaches p=fjLs when x is large and positive, or p— — fis when x is large and negative. 

This lowest point of exp{—px)y[pyUxk) against real values of p constitutes the one and only 
saddle point along the steepest descent path when ex^{—px)y{p,Ux, X) is integrated with respect 
to the complex variable p (see fig. 3). 

The minimum and the saddle are sharp when \x\ is large. Under these conditions, a knowl- 
edge of y{p,Ux, X) over a narrow range of real values of p suflSces to evaluate the integral (7) 
and thus to calculate the total distribution of X-rays far from the source. The main object 
of this paper is just to study the trend of y for real values of p near /Zs (or — Ms). Section 9 
deals with the actual evaluation of the integral (7) . 

For smaller values of the depth of penetration x, the saddle point lies nearer to p=0. 
The direct calculation of Y(x,Ux,\) at moderate depths also provides convenient expansions 
of?/ in powers of p, or of related variables like p/(fJLs—p) and p^/ifh—p^) [13]. This infor- 
mation may serve to interpolate y between the proximity of p=0 and of p= dzMs- 

5. Special Analytical Solutions — Straight-ahead Approximation 

In order to study the trend of y for real values of p near fig, we seek approximate analytical 
solutions of the transform equation (8) valid for values of X'^Xs and of Ux^^l. This search 
will be conducted first in the straight-ahead approximation even though this approximation is 
inadequate. 
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The straight-ahead approximation consists of setting Ux=l on the left side of (5) or (8). 
This procedure is equivalent to the replacement of the depth of penetration a^ in (5) with the 
path length traveled by the photons, irrespective of the obliquity of penetration [18]. The 
variables Ux and u that still remain in the eq (8) can then be eliminated by integrating over u. 
Setting 

I duy(p,Ux,\) = Hp,\), (9) 

we obtain the straight-ahead approximation equation 

[fi{\)-pmp,\)= rHy,\mpX)dv+8i\-\o)i, (10) 

where 1= f{u^du. ^ 

J47r 

Because we are interested in an approximate solution valid over a narrow range of variation 
of X, we disregard the variations of the differential scattering coefficient k over an interval 
X''^X and take 

A:(X',X)-A»(X,X)=(3/4)/ir, = 0, (11) 

in agreement with the definition of k (see section 3). 

Before entering this approximation into (10), we must distinguish two different situations, as 
in section 2.3. (a) If the primary radiation is the most penetrating one, that is, if iUs=M(^o)=Mo, 
the range of approximate integration over X is right next to the source wavelength Xq. 
(b) If the most penetrating radiation is that of wavelength \m, for which ^ has its absolute mini- 
mum value iJL„i=fi{\m)j that is, if }^o<C^m, the range of approximate integration cannot run all 
the way from Xq to X^. In this event, we do not attempt to solve the inhomogeneous equa- 
tion (10), but only the corresponding homogeneous equation, in the range of X^X^„. 
Therefore, we seek the solution of 



[^(\)-pM\p)=C rH\\p)dy+1 6(X-Xo) 



for 

X-Xo>X,, (12a) 

and of [fi(\)-p]^(Xp) = Cr^{r,p)dy 

for 

X~X,,>Xo. (12b) 

These equations reduce to differential form, if a derivative is taken with respect to X, and 
their solutions are 

HKp)=I(^o-p)-'C [^(\)-p]-' exp [cf^dyMV)-p]^ (13a) 

H\p) = C[fx{\)-p]-' exp |cJ'rfX7[M(X0-i>]}- (13b) 

These solutions are derived as special cases of a more general procedure in [18 section lib]. 

We are interested in the behavior of these solutions when p approaches the singular point 
Us, that is, juoin (13a) and /x^in (13b). 

In case (a) the first derivative of m(XO ^t the singular point /xo, (^m/^XOo^mo, is different 
from zero. Therefore, the integral in the exponent of (13a) has a logarithmic singularity. 
Incase (b) the first derivative of /x(X') vanishes at iu(XO = Mm, X' = X^. Therefore the integral 
in the exponent of (13b) over an interval of X'-^X^ is of the arctangent type and its singularity 
has the form {yLm—p)'^'^. If m(^0 is regarded as linear in case (a) and parabolic in case (b), 
the solutions (12) take the characteristic singular forms [3]: 
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#(M~I(mo-7'^-'<^[m(X) -;;]-• exp {-2 log ^^^^^j=IC'[M(X)-;'r^"-'[/x,-y>]-^''^o-- (Ha) 

(Mo Mo~P ) 

H\p)-C [uW-p]-' exp {27rC/[2M,n(M. -i>)l'^^'}, (14b) 

where ilm = {cP^Ji/d\%m' 

Even though the sohition of the transform equation (10) for vahies of X>>X.<f cannot be 
obtained in this simpk^ manner, its dependence on p, as p approaches ijls, is described by (14) 
for the following reason. As soon as X is sufficiently larger than X^ so that /x(X)— ^>>/Zs— i>, 
p may be replaced in the equation by /i„,. This implies that the further variations of <l>(X,p), 
as X keeps increasing, are independent of ^, that is, that a state of limited equilibrium, within 
the meaning of section 2.2, has been achieved. As p approaches Hs, X need depart from X^. less 
and less, before the condition /x(X)— p>>Ms— i>is adequately fulfilled. Therefore the approxi- 
mations made in the derivation become increasingly good. 

The inversion of Laplace transforms with the singularities (14a) and (14b) will be con- 
sidered in section 9. However, we anticipate here that the inverse transforms behave, respec- 
tively, as (2) and (3) for large depths of penetration x. 

For the sake of completeness, we mention yet anotliei' situation, which combiners the 
features of (a) and (b), namely, the case where Xo<X<X„„X,s = X. Here the pertinent form of 
the equation is (12b) and its solution (13b). However, the derivative (d/ji/dX),, at X^X., is 
difl^erent from zero. Hence, the integral in the* exponent of (13b) has a logarithmic singularity 
and we find, instead of (14b), 

</)(X,;:>)-C[m(X)-2>]^/A«-i. (14c) 

Various attempts were made in the past to develop methods of succ(^ssiv(^ approximations, 
starting from the approximate solutions (13) and (14), to solve the straight-ahead equation (9) 
for all values of X. Direct numerical integration eventually proved superior [6], even for values 
oi p rather near ^m despite the fact that the rapid rise of y along X, for X near X.„ requires small 
steps of numerical integration. 

6. Possible Treatments of the Obliquity of Penetration 

As pointed out in section 2.3, the deep penetration of X-rays is influenced to a substantial 
extent by the small obliquity of the most penetrating components. In order to evaluate this 
effect, one must, in principle, work out the directional distribution of the X-rays of each wave- 
length, that is, solve the complete two-variable equation (8). In practice, one may want to 
learn about the distribution-in-angle only the least amount necessary to evaluate the eftect on 
the trend of penetration. 

In order to analyze this situation, one may expand the distribution-in-angle of y(p,Ux,}^) 
into Legendre polynominals Pi(Ux), 

y(p,Ux,\)=J2ii2l+l)(47r)-hji{p,\)Pi{Ux). (15) 

with 

yi(Pj\)=j^^duPi(u,)y{p,u,,\). (16) 

Equation (6) multiplied by each of the Pi^s and integrated, reduces to the system of equations 

M(X)2/o(2>,X)-Mi(:P>) = £/"(X',X)2/o(2^,XOrfV+5(X-Xo)£/(^.)rfM 
(17) 

M(X)2/K2>,X)-i>[(/+l)i/z+i(p,X)+/2/,_i(p,X)]/(2/+l) 

-J^'^/:(X',X)P,(l-X+V)t/,(;>,V)rfX'+6(X-Xo)J^^P,(uJ/(i^,)rfM. 
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The first equation of this system has the same form as the straight-ahead equation (10), 
except for the fact that t/i appears in the second term on the left side instead of y^. If there 
were a way of guessing the ratio 

]/i(i>,X)/yo(p,X)== r u^y{:PyU^yX)dul\ y(p,u,,\)du=u,(p,\), (18) 

yi could be replaced with Ux{p,^)yo(p,^) in the first eq. (17), which could then be solved as a 
single-variable-straight-ahead equation. Thus a knowledge of Ux(p,^), or of an equivalent 
parameter, embodies all the information about the directional distribution that is immediately 
relevant to the penetration problem. 

Instead of Ux one may want to use the equivalent parameter [6] 



gi 



(p,\) = l—uj^p,\)= I (l—u,)y(p,u^,\)dul\ y(p,u^)dUy (19) 



which represents the first moment of the angular distribution. If we set yi=--{l-gi)yo in (17),. 
the first equation becomes 

[fjiW-p+pgi(^,\)]yo{p^)= rk(W,X)yo{p,\')dy+d(\-\o) f JMdu. (20) 

Thus, gi(p,X) is seen to measure the departure of the eq (20) from the straight-ahead eq (10). 

In the critical range of wavelengths, X^^X^jthe solution of (20) may be expressed in a manner 
similar, for example, to (13b), that is, 

^o(^,X)-C[M(X)-:p+p^i(^,X)]-iexp|c'JW/[M(X0-2>+m(pA0]| 

^+vgi{\v)IW{y^)-p\^\j i+?>^i(p,xO/[m(xO-p] m(xO-p) ^ ^ 

The last expression differs from (13b) in that C is now divided by l-{-pgi{\,p)j[ii{\) — p\. We 
shall see in the next section that this factor tends to a constant limit as p approaches the singular 
point Ms- 

If the ratio Uj^—yjy^ is assumed to be known, y^ can be calculated, and thereby yi becomes 
also known. The second equation (17) yields then 2/2 directly, the third yields yz, and so on, 
that is, the whole system of equations unravels automatically. It was pointed out by Wick [1], 
and also by Waller [19], that the set of values thus obtained for the ?/z, diverges, in general , 
unless the correct value of Ux was chosen at the start. Thus the solution of the system (17) 
takes the character of an eigenvalue problem. (This property is common to all systems of 
which the first n equations contain n-{-l unknowns). 

A simple trial and error procedure on the choice of Ux=yi/yQ in (17) is impractical in the 
critical range of X where the directional distribution is sharply peaked forward and the successive 
coefficients yi,yi+i . . . are nearly equal, so that the trend of the sequences yi evolves exceed- 
ingly slowly. In other words, the angular distribution cannot be characterized in this range 
by just a few Legendre coefficients. Therefore Wick [1] replaced the sequence of the ?//s 
with a continuous function and the expression [(l+l)yiJ^i+lyi-i]/{2l-{-l)—yi in (17) with the 
second derivative of the function. In this manner the problem of determining the value of 
"^x^yilyo was reduced to a problem of differential equations. An approximate determination 
of the eigenvalue by variational methods offers little difficulty and was carried out by Wick for 
a constant m (X) and for a linear /x(X) [1]. 

Whereas variational methods are well suited for deriving rapidly an approximate eigen- 
value, they tend to yield poorer approximations to the eigenfunctions and generally do not 
appear convenient as a basis for higher approximations. In fact, Wick effectively characterizes 
the angular distribution of radiation by a single adjustable parameter, in the first approxima- 
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tion. A substantially better approximation is required to take into account more gejieral 
variations of m (X) as well as variations of the kernel k (X', X). Extensions of Wick^s procedure 
in this direction have been developed but they appear cumbersome. 

When the directional distribution is sharply peaked, as it is in the critical range, X'^Xs, it 
can be characterized conveniently by its set of moments, which in this case tend to converge 
far more quickly than the set of Legendre coeffcients yi. The moments, mo=2/o) '^i=9iyo= 
Vi—yoj . . . 7nn=\ (1— Ua;)"2/(Pj^a: A) ^w are linear combinations of the Legendre coefficients. 

Determinations of ^1 by trial and error procedures based on calculations of the m„/s should be 
more practical than procedures based on the y/s,. 

Spencer eventually developed an efficient iteration procedure that permits to determine 
^1 by trial and error, utilizing only a very few^ — as few as two — of the equations for the 
moments [6]. This method has superseded all previously attempted approaches. 

One point, which emerged in the course of earlier attempts, might be worth mentioning 
because of its possible more general interest. The testing of trial values of ^i(p,X) becomes 
particularly easy when the obliquity of penetration is less effective than the variations of 
/z(X) in shaping the X-ra}^ spectrum, that is, according to section 2.3, when d log /z/(iX<^l. 
In this event, an incorrect initial assumption regarding gi, not only causes the sequence of 
moments mi,m2 . . . m^ . . . eventually to diverge, but quickly to assume a grossly erratic 
trend. This critical instability appears to be a general property of systems that may be 
called ''strongly nonself-adjoint." It is discussed in appendix B. 

7. Effects of Small Deflections on Special Analytical Solutions 

As mentioned in section 6, Wick did the basic work on the effect of small deflections on 
penetration [1], for two typical situations, namel}^, constant ^(X) and linear ju(X). Either of 
these assumptions about m, together with the assumption (11) of a constant kernel ^(X',X), 
makes it possible to separate the direction variable u^ from the energy variable X in the transform 
equation (8). 

Specifically, the homogeneous part of the equation thus simplified has solutions that are 
products of a function of X and of a function of (l—Ux)p/[iJi(\) — p]. This to say that the direc- 
tional distribution maintains a constant shape as X varies, but it contracts in width, as fx{\)—p 
grows smaller, in proportion to [m(^)— ?>]/?>• As a result, the factor i+p^i(X,p)/[M(X)— p] 
in (21) remains independent of X and p. Therefore, the dependence of the photon 
number on X and p, irrespective of direction, which is given by (21), has the same analytical 
Jorm as the straight-ahead approximation (13b), except that C is replaced with the constant 
C/{l-i-pgi/(fjL—p)}. There is a whole set of solutions of the homogeneous equation, with 
difl'erent values of pgi/(fJi—p). 

Wick superposes solutions of this set so as to construct a solution of the inhomogeneous 
equation with the 5(X— Xq) source. The component solution with the smallest value of 
pgi/{fji—p) predominates over the other components as p approaches /x^ and as X increases 
through the critical range X-^-Xs. Therefore, the singularity of the Laplace transform is char- 
acterized by the variation of ^(X) — p in (21) (or in (13b)) as p approaches fis and by the smallest 
value of pgi/{fJL—p). 

The application of this procedure to the X-ray problem is straightforward when the pri- 
mary radiation is the most penetrating one (Xo>X^i, /Xs==Mo) [4], since in the very proximity of 
Xo the variations of /z(X) may well be regarded as linear, m(^) '^Mo+Mo(>^— ^o)j and those of 
k(y,\) may be disregarded. 

The applicability is less clear when the primary radiation is not the most penetrating one 
(\o<C^m), that is, when fi{\) goes through a minimum in the critical range. However, one 
may surmise that the problem at the minimum (X^X,^) can be treated by assuming the same 
value of pgi/(iJL—p) as though ju(X) were constant. This assumption, which has not been dis- 
cussed previously, implies that the singularity of the transform is given by the straight-ahead 
formula (14b) with Creplaced with the vahie of C/{l + pgi/{n—p)} corresponding to a constant jjl, 
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We present here first a qualitative argument and then a more detailed treatment which 
confirm the correctness of the surmise. 

The trend of the solution of the transform equation in the critical range depends on the 
variations of the attenuation coeflBLcient /x(X) and on the efl^ect of small deflections. We want 
to show that the deflection effects predominate over the variations of m(^)j in the limit as p 
approahces ^im^ so that the deflection effects may be treated as though /i(X) had the constant 
value ixm- 

Obliquity of penetration effectively impedes the build-up of photons if the value of 
iJL—pUx in (7) becomes substantially larger than the corresponding ^ ^straight-ahead'^ value 
iJL—p. Therefore, the main build-up is confined to photons for which l—Ux'^{iJL—p)/p. 
Compton scattering of photons within this range of directions corresponds to wavelength 
shifts 6X< 1 —Uj,< {iJL—p)/p. 

On the other hand, variations of the absorption coefficient that affect the asymptotic 
behavior must be of the order A^i^firn—p- These variations may be described, in the neighbor- 
hood in the minimum, as Ajli'^Jju^AX^ (where iJLrn={d^t^ld\^)m)' Therefore, the wavelength 
change AX required to bring about a substantial change of jjl is of the order of magnitude 

Hence, in the limit for iXm—p—^O, this shift AX becomes proportionately much larger than 
the shift d\ required to bring about substantial effects of obliquity. 

Proceeding beyond this qualitative argument, a complete treatment of the problem would 
consist of solving the transform equation (8) by a procedure of successive approximations, 
whose first step should yield the Wick- type results. Efforts in this direction have yielded 
rather cumbersome developments. In practice such a complete treatment has been made 
unnecessary by the development of the semiasymptotic numerical method [6] (see section 2.4), 
which requires, as an analytical complement, only a knowledge of the singularity at p=fJis, 
Ux=l. One possible approach to the remaining problem would be to seek a special solution of 
the eq (8) for the case when the primary radiation is not the most penetrating one, that is, 
when \o^^m, i^s=i^m, since Wick has already solved the problem for Xo>X;,i, Ms=Mo. Whereas 
Wick assumed a linear variation of /i, near X = Xo, we should assume here a parabolic variation 
near X^X^^^, that is, M(^)'^Mw+2Mm(X— X^^)^. Instead of following this approach, we shall take 
a somewhat more general one, namely, to seek special solutions of the equation valid ^^in the 
asymptotic limit'' 

[^W-pVp^O, (22) 

that is, for X^^X^, Ux^^l. This will be done without utilizing any initial approximation about 
)u(X) or k{y,\) but the approximations considered above will result automatically from the 
consistent application of the limit process. 

The first step of our treatment consists in choosing instead of u^ a new direction variable 
that measures the effect of obliquity on penetration. Equation 21 shows that this effect is 
proportional to the mean value of p{l—Ux)/[fJi{\)—p]. Accordingly, we replace the unit space 
vector u={'d-,<p) with a homolographic projection variable that is represented in plane polar 
coordinates by the vector 

»=([2p(l-tt,)/(M-2^)]"^^) ~([W(M-i?)r'^«?,^). (23) 

(The approximate equality corresponds to the small angle approximation ??<^<^1, which Ave 
need not utilize at this point.) The inverse transformation of (23) is 

Ux=l — W[yiW—p\lp, u(v,\) = (sirQ.osUx,(p). (24) 

It is also convenient, though not essential, to renormalize the dependent variable y{p,Ux,\) 
by multiplication with iji(\)—p, that is, to take 

w(v,\) = [m(X) -p]y(:PjUx,\). (25) 
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The dependence of iv on p will not be indicated explicitly. • 

The transform equation (8) takes now the form 

{i+^v')w(v,\) = {27rp)-' i (/yk(y,\) fdv'd(\-U'U-\+y)w{v',y)+id(\-\,)j{u,) 

=={2Tp)-'jdv'k(\-A[\v,v^\)w{v',\-\[\v,v'])+ld{\^^^^ (26) 

where A stands for 1 — m-m'. Since u is given by (24) as a function of v and X, and similarly 
u' as a function of r' and X' = X — A, A(\,v,v^) is defined as the root of the ec{uatioii 

A=l-w(r,X)-M'(r',X-A). (27) 

We are now interested in the solutions of (26) in the critical region X^X^, /x(X)^/x.s and in 
the asymptotic limit [fi{'\)—p]lp'^0. Following Wick, w^e seek here solutions of the homo- 
geneous equation corresponding to (26), that is, of 

(l+^r')w{v,\) = {27rp)-'j(Iv'k{\-\,\)w(v',\-\), (28) 

and leave it for the next section to superpose solutions of (28) in such a maimer as to fulfil 
the inhomogeneous equation (26) at the source where X=Xo. 

Our knowledge of the straight-ahead solution and the discussion of section 6 on the effect 

of small deflection suggest that we split off from w(v, X) a factor exp IC I ^/X7[m(^0~":P] (' 

where C remains to be determined, with the intention of showing that this factor contains the 
wdiole singularity of w m the asymptotic limit. Thus we set 

w(v,X) = x(v^}^ exp \cj\iyMy)-p]] (29) 

y(p,u,,\) = M\)-p]-^x{v^) exp \cj\iyMy)-p]^' (29') 

Equation (28) becomes 

(l + h^')x(v,\) = (27rpr'jdv'k(\-A,\)exp\^-CJ^_J^^ (30) 

and we shall show that in the asymptotic limit A vanishes, so that k(\—A,\) equals C, and the 
exponential becomes independent of X. As a result x becomes, in this limit, a nonsingular 
function of v only, so that the exponential in (29) does indeed represent the singular behavior 
of w. The values of C will be fixed by the condition that x (^) converge for large v. 

To carry out this program, we begin by studying the equation (27) which determines A, 
If we substitute u and a' from (24), (27) becomes 

Terms of higher order in [ijlO^)—p]/p will be disregarded henceforth in the asymptotic limit. 
Such terms should be treated by a power expansion in case one intended to proceed to highei 
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approximations. This procedure has the same effect upon (31) as the ^'small angle approxima- 
tion'^ indicated in (23). Equation (31) reduces then to 

A=l ^^(^^ ir-^+^^^^7=^ -'1-r^^%^^T"2 v.v'\. (32) 

2 ^ (L f^W-v J L f^W-p J ) 

Since this equation indicates that A tends to vanish in proportion to [fJLO^)—p]/p, in the 
asymptotic limit, we replace A with the variable 

^=ApM\)-pl (33) 

which we surmise to remain finite in the asymptotic limit. The dependence of ai(X— A) on A 
will also be indicated through a parameter which tends to a finite limit, namely 

r7^[M(X)-M(X-A)]/>A==[M(X)-M(X-?[M(X)-i^]/:p)M/x(X)-;>l (34) 

M(X-A) = M(X)-r7pA. (35) 

With these substitutions, (32) reduces to the quadratic form 

(2+v'\y^'~2[2iv'+v'^) + (y' + vV'-2V'V'')rj]^+[(v'+v''^^^ (36) 

This equation depends on ^ not only explicitly but also through 77. However, if the right side 
of (34) is expanded into powers of [iJ'(X)—p]/p and only the first term is retained, one finds that 

V^.{d^/d\)lp = ijiW/p (37) 

independently of ^. Therefore, f is given by (36) as an algebraic expression ^(v-v^ii/p) which is 
finite and independent of M^)—p]/p in the asymptotic limit. Owing to (33), A vanishes then 
in the asymptotic limit, as we had surmised. 

Therefore, in this limit, k(\— A,\) =k {\,\) = C, the exponent of (30) becomes 

C r dyMV)-p] = C r dyMV)-p] = {Clix) in [l-^i^M. (38) 

Jx-A Jx-^[m(x)-p]/p 

and the equation (30) reduces to 

(1 + W)x(v^) = iCI2Tp) jdv' [1 - ? (v,v',ij./p)ii/pf':'x(v',\) . (39) 

The kernel of this integral equation depends on X through jjl. Within the narrow critical 
range of photon energies, X-^X^, in which we seek an anah^ical solution, fi (X) can be expressed 
as a power series 

m(X) = A..+ m.(X-X,0 + 1/2m;(X-X,)'+. . . (40) 

Accordingly, one may seek the solution of (39) by expanding both the kernel [l — ^(v,v\ ii/p) ixlp\CJii 
and the unknown function x(^A) into powers of (X— X^). Since the critical range around Xs 
becomes infinitely narrow in the asymptotic limit, the singularity of y(p,Ux,\) (at p=iJLs,Ux=l) 
depends only on the zero-order term of the expansion into (X—X^)"^. Therefore, for the purpose 
of studying this singularity, we may set 

x(v,\)-^x{^y\) = Xo{v)y (41) 

where xo(^) fulfils the eq (39) at X = Xs, and with p=fJLs, that is, 

(l + iv')xo(v) = {C/2Tfji,)jdv^[l~^{v,v\ijiJiJis)ijis!^sr^'^ (42) 
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In this manner we have reduced the problem to the solution of an equation in tlie sino:I(^ vari- 
able V, that is, we have achieved the same separation of the energy and direction variables 
(X and v) as was achieved by Wick. The only difference lies in having achieved tlie separation 
by a limiting procedure rather than by initial assumptions regarding the variation of /z(X) and 
the constancy of k(\,\). Notice that (42) holds whether X^ coincides with the primary wave- 
length Xo (that is, for Xo>X;,J, with the longest wavelength X (that is for Xo<X<X;,,) or with the 
wavelength of minimum absorption X,„(that is for Xo<CXw<C^)- Iii the last event iis = jx^=0 and,, 
according to (36) and (37), the kernel of (42) reduces to 

{[l-f(l;,l?^A.7M.)AJMJ'^^2^•]M.=o==exp[-(^/MJ?(^,l?^0)] = exp{-K^^^ (43) 

Equation (42) is equivalent to the Wick eigenvalue equations and has finite solutions xo(^) 
only for special values oiC, As in the standard eigenvalue problems with a single independent 
variable, the successive eigenfunctions Xo^^(^),Xo^^(^)j . • • Xo^^(^) • • • have 0,1 . . . n . . . 
nodes. The corresponding eigenvalues C^^\ C^^\ . . . C^"^^ . , . depend monotonically on n. 
The eigenvalue C -^^ is the largest one, algebraically. Each of the eigenvalues C^'*^ characterizes a 
type of singularity of y{j),Ux,\)- As mentioned before, these singularities coincide with those of 
the straight-ahead solutions (14) in section 5 except for the replacement of the constant C with 
(7. The solution of the inhomogenous equation (8) or (26) will be constructed, in the next 
section, as a superposition of components, each of which is a solution of the homogeneous equa- 
tion, with different values of C. We anticipate here that the component with the largest value 
of C in (29') is clearly the one with the sharpest singularity and therefore the one which predom- 
inates and which alone matters in the asymptotic limit. 

The equivalence of (42) with the Wick equations is immediately apparent when \s = ^mr 
l^s=f^m, f^s=^, that is, when the kernel of (42) takes the form (43), so that 

^l+lv''^Xo{v) = [C/2n^jfdv'exp^-^^ (44) 

This equation coincides with (W53) [1], as shown by substituting 

v=^{M!ky^'s, C/fx^n.=2k/M, Cly.m={M+\yi2M\ (45) 

When Ms 7^0 no better procedure is seen, to show the eigenvalence with the Wick equations, than 
to backtrack from (42) to an equivalent homogeneous equation in the variables X and Ux and ta 
handle this problem by Laplace and Fourier transforms in the manner of Wick. One may then 
identify xo(^) with a transform of solutions of Wick equations according to the formulae 

Xo(z;) - (/i,/27riM.) {'^ dt exptt'-'^^^sjmiisjfXsY^H) 

J —ico 

= (!is/8irH) C" dtjdqexp (t-if'^q-v) <'-^'^0([M./2Ms]"^g) 

= {its/8irH) f" dt (dq exp {t-it'"q-v+qyu)f-^''''.U(M2-^i,yJ'q). (46) 

Tlic functions/, $ and U obey the equations 

^-{i-^lLMd/du)+\^,u'y{u) = {C/in)jdu' exp (^ I «--«?)/(«'), (47) 

^^nA(dldcy+<r-'(d/da)]-^i,Md!da) + Cex^ (- c')^4.{c) =U4>ia) , (48) 

^jf,A{dld<Tr+cr-'{d!da)]+^^,+C exp i-<T')-~i^y/^.s} U{a)=CU{c)- (49) 
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All of these equations obtain when /Xs=0 as well as when Ms?^0, aiul for Ms<CO (that is, for Xo< 
Xs<Xm) as well as for /X5>0. In particular, (47) coincides with (44), for ais=0, if u is replaced 
with (20 Ifx^y^'v. 

The self-adjoint differential equation (49) has the form of the Sclu^oedinger equation for 
the motion of a particle oscillating radially in a cylindrically symmetrical potential. It is 
identical watli (Wl05) [1] as sliown by performing the substitutions 

(7/M.= (M+l)V2M^ iiJfi, = 2y/M, C/i^s=2-p. (50) 

It is also identical with eq (V4) [4] as shown by replacing a'^ with 0-72 and C with —QUxs- 

Wick has estimated the largest eigenvalue of his Schroedinger-type equations by a vari- 
ational method. When /is=0 he finds 

c^{C'''-h^j''f=c[i-h{tj.jcy''Y. (51) 

This eigenvalue reduces to C as ^JC approaches zero. The ratio t7//x,^ represents the differ- 
ential cross section for Compton scattering, for X' = X, expressed in units of the minimum value 
of the total cross section. When /is^O the largest approximate eigenvalue is 

^^M.[(C/M.)r(r+l)-i + iA./M.-i(A.V.)-V-^ (52) 

where the parameter r indicates the solution of the algebraic equation 

4(6^/m.) (r+ 1)-H (m./m.)V-^== 1. (53) 

When the logarithmic derivative |/xs/msI is large, one may solve (49) by a perturbation 
method starting from the solutions of the quantum mechanical harmonic oscillator problem 
and treating e = Hs/ij.s as a small quantity. The result, given in (V6) [4], is 

C=--C-hiji(a-l) + Ceya(a + e) + Ch'(2a + ey/iji,a'{a + ey+ . . ., (54) 

with 

« = [1+4(C/A.)e]". (55) 

In practice idslids is seldom very small, and the range of application of the simple perturbation 
method is limited. The investigation of eigenvalues proceeds more effectively by the simple 
variational method of Wick or by the methods that were developed by Spencer [6] for the numer- 
ical solution of the semiasymptotic problem. Appendix C contains a description of these latter 
methods and tables of eigenvalues. 

8. Effect of Obliquity of the Primary Badiation 

We resume here the discussion of section 2.6 regarding the influence of the source direction 
on the penetration law. 

In the last section, w^e have studied the singularit}^ of special solutions of the transform 
equation at p = iis- These special solutions obey a homogeneous transform equation shorn of 
any source term. It was understood that the solution y{])^Ux,\) of the complete, inhomogeneous, 
transport equation (8), in the critical range X'^X^ and in the asymptotic limit, should be con- 
structed as a superposition of special solutions (29') 

7/(^,7/„X) = SA[/x(X)-p]-^X^"U^,X) exp{^^^)J'(/V/[M(V)-p]} 

^Z,aM^)-V]-'x^\i^ exp|c(«>J'rfV/[M(V)-i>l|- (56) 
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Here xd^'^C^O a,nd C^''^ represent, respectively, the 77--tli eigenfunction and eigenvalue of (42). 
The function x^''^^',^) in the middle part of (56) represents the complete expansion 
X^'*^ = 2;.Xr''^(y)(X— X,)^ which begins v^itli xd'"'* and can be carried to highei' ordei* of ;ij)pi()\iina- 
tion if necessary. 

As^ approaches ju,, the factors exp< C^""^ I d\' l[^{y) — j)\ > become infinitely large, ])rovi(hHl 

C^^'^^O. Wick pointed out that only one of the eigenvalues (7^"\ namely C^^\ may happen to 
be positive. At any rate the term n=0 of the sum (56), with the algebraically largest eigen- 
value C^^\ will increase fastest as j) approaches ai^ and will predominate over all others in the 
asymptotic limit. Therefore, the singularity of y{p,Ux,\) depends primarily on the behavior of 

exp < C^^^ ^X7[m(X0~p] [ > ^s anticipated in section 7. The behavior of this factor is described 

by (14) in section 5. 

In addition, one must consider the factor [m(X)— ^p]"^ if ix{\)^}is and especially the depend- 
ence, if any, of the coefficient a^ on j). The coefficients an of the superposition (56) are deter- 
mined by the inhomogenity of the transform eqn (8), that is by the characteristics of the radia- 
tion source. Therefore, these characteristics may influence the singularity of y{p,Uj:X), and 
hence the penetration law at great deptlis, through the dependence of a^ on p. 

It has been pointed out in section 2.6 that the obliquity of the primary radialion has much 
less influence on the penetration for high source energies, Xo<C<C^/n than for Xo>X.^j. This im- 
portant difference shows up clearly in the process of determining the coefRcients an of the ex- 
pansion (56). 

If the source energy is considerably larger than the energy of maximum penetration, one 
will probably have to resort to numerical integration of the transform eq (8) from the source 
wavelength Xo at least up to some wavelength X/ somewhat lower than X,,^. If \<C^\n, and one 
is interested in the distribution of photons of energy higher than the energy of maximum pene- 
tration, then Xs equals the wavelength of these photons and one must choose X/<Xs. Whether 
Xs<C^m 01' Xs=X,„, we define X^ so that \<i\f<\s^^m- A superposition of analytical solutions 
valid in the critical range X'^X., should then be fitted to the numerical solution at X/. Under 
asymptotic conditions, (m^— |>>)/;^<<1, the fitting-wavelength X/ can be so chosen that m(X/) — 
jLts^^ju^— 2>, that is so that ix{\)—p can be safely replaced with fi(\) — fjLs for X<X/. When this 
is so, the photon distribution at X/ varies no longer as p a])])roaches jLt^ still furthei*. Similarly 

the analytical solutions x^"K^'A) exp<(7^''^ ^^^7[m(^0~P] [? extended from X^X^ to X/, no 

longer depend on p in the proximity of X/, where p may be effectively replaced with ijl,. 

Under these conditions, the fitting of the superposition of analytical solutions to the actual 
photon distribution at X/, that is, the determination of the coefRcients a« of (56), gives a result 
independent of j^- Any change of obliquity of the source radiation causes a change of the 
directional distribution of the photons at X^^ and, therefore, a change of the coefRcients an. The 
photon distribution through the critical range will vary with p in a manner which depends on 
the intensity ratio of the component solutions anXn(i',^) exp { . . . } of (56) and, tluMu^fore, on 
the obliquity of the source radiation. 

However, the intensity of the lower energy radiation depends only on the n=0 term of (56), 
once p is sufficiently near to fim, and therefore it will vary with p only through exp 

\C^^^ J dyi[ix(y)—p]\, independently of the ratios among the a^'s. That is, the source obliq- 
uity influences the intensity of the lower energy radiation through the same factor ^o independ- 
ently of depth of penetration, provided this depth is adequate to insure fully asymptotic 
conditions. 

A quite different situation prevails when Xo>X;„, since the coefficients of the analytical 
solutions in (56) must be fitted to the directional distribution of the source radiation right at 
Xo= Xs, in the middle of the critical range. Here, the dependence of the anS on p has been worked 
out by Wick for an isotropic source [1], and in reference |41 for a monodirectional source. 
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General formulae equivalent to a determination of the anS will be given here below, but they 
are of very limited usefulness because of their complication and, especially, because they imply 
a knowledge of the directional eigenfunctions Xo"^('^)- A practical solution of the effect of 
source obliquity under these conditions will probably require much additional work. 

An analytical solution for Xo^X;;^, equivalent to (56) and fitted to the directional distri^ 
bution of the X-ray source, can be obtained performing a Laplace transformation of the energjr 
variable X, according to Wick [1]. This is a solution of the transport equation (8) schematized 
in the following manner: 

(1) jLi(X) is taken as mo+Mo(X— Xq), 

(2) k{y,\) is replaced with C, 

(3) u= {^,(p) is replaced with a ^'small-angle approximation'' vector s= (s,(p), with <s < oo ^ 
such that Ux may be replaced with 1 — Js^, and wu^ with 1 — | |s— s'|^. For small angles, one> 
sets S'^^. The solution of the equation is then written as 

1 r^°° 1 r - 

ZttI J -im 47r J 

X^JJn{k)^L- rdle-^' fl+^l''"'"" D^ (r, ri + ^> (57> 

Mo — pjo L Mo — pj V L Mo — Pj/ 

Here Un and C^^^ are the eigenfunctions and eigenvalues of (49) and 

is, in essence, a coefficient of the expansion of the directional distribution of the source — in- 
dicated hjj{ux) in (8) — in terms of the eigenfunctions of (47). Some simplification is attained 
by considering only the X-ray flux integrated over all directions, since in this case (57) yields. 

(59> 

Much further simplification obtains for two special types of source, namely, (a) mono- 
directional perpendicular to the source plane,/(l — |«^) = 5(s), and (b) isotropic, or practically 
isotropic, /(I — |s^)'^7=const., since (58) reduces, respectively, to 

D,(a)=^ r kdke-^ok^l^Wn{k), (60a> 

^T^ Jo 

and _ 

D^(a)=iUniO)f. (60b> 

Whether Dn is constant or proportional to its argument, the two integrals in (59) may be^ 
treated as a Laplace transform followed by its inverse, so that (59) reduces to 



and 



J,J,!,(p,8,X) = 2.1[R(0)r/^[l+S^^J'"'''^-', (61b> 



respectively, that is, to algebraic functions of mo— P- (The p in the exponential of (61a) may ba^ 
safely replaced with /jlq) . 

For other types of directional distribution, the dependence of (57) or (59) on fiQ—p is quite: 
complicated and even difficult to survey qualitatively, as had been anticipated. 
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9. Inversion of the Laplace Transform 

Following the discussion of the singularity of y(p,Ux,\) at p==Ms ^^ sections 5 and 7, we 
shall now translate each type of singularity of the Laplace transform into a corresponding 
distribution-in-depth of the X-rays for large values of the depth x within the medium. The 
translation is done by means of the inverse transform formula (7) 

I poo ^ pc» 

Y(x,Uj,,\)=:z—. exip{—px)y(p,u^,\)dp = exp{—fisx)w—' exp [{^s—p)x]y{p,u,,\)dp. (62) 

ZttI J —ioD ^7^^J-^oo 

Having split off the factor exp( — ju^x) in the last expression, the remaining integral represents 
the buildup factor B{x,Ux,\) (to within a constant factor equal to the source intensity). Hence • 
forth, we shall be interested in the dependence of J5 on a:: for any given obliquity Ux and wave- 
length X and, therefore, we shall omit explicit reference to u^ and X. Similarly, it will be con- 
venient to indicate y{p,Ux,\) simply as yifj^s—p), so that (62) takes the form 



1 r^°° 

^TTl J —IOC 



B(,r)=^^. dp exp [(iJis-p)'r]y(fx,-p). (63) 

The preceding sections have called attention to two types of singularities of y(iJ's—p)i 
namely: (a) an inverse power law, which we indicate as 

y(n,-p)cc(iJLs-p)-^-\ (64a) 

and (b) an exponential function of the inverse square root, which we indicate as 

y(iJLs-p) oc exp [DfXs/(iJis-p)Y^K (64b) 

Equation (64a) describes the singularity when Xq^X^^^ and the directional distribution of the 
source is either perpendicular to the source plane or effectively isotropic, as shown, respectively, 
by (61a) and (61b). In addition, the same type of singularity arises for \o<C\n in the calcula- 
tion of the buildup factor for X^X^, that is, when Xo<CXs<CXw, in which case the straight-ahead 
result (14c) applies, except for the replacement of O with C^^\ Henceforth C^^^ will be in- 
dicated briefly as C. The value of K in (64a) is given by 

K=C/\ix,\ (65) 

for \s<C^m, according to (14c), and for }^s^\n with an isotropic source, according to (61b), 
whereas 

K=C/iis+l (66) 

for ^y^^m with a monodirectional source, owing to the derivative in (61a). Equation (64b) 
describes the singularity when X^^^X^ and X>X^, which implies Xo<X;^, in which case (14b) 
applies, except for the replacement of C with C, and 

D=2^'C'lfjiUd'n/d\')m. (67) 

It may be added that the build-up of the radiation of maximum penetration X==X;^, for 
which \s=\^ also,* involves an integration of the arctangent type of integral in (13b) over a 
range of 7r/2 of the arctangent, instead of t. Moreover, in this case, the factor [m(X)— 7>]~^ in 
(14b) represents (ms— 1^)"^ As a result, (64b) must be replaced with 

y(iji,—p)oc{fjL,—p)-^ exp [FiJLsK^s—p)V^\ (64c) 

where 

F=Tr'Cy2iiJ4\'^ld\')m=Dl4. (68) 
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Notice that both K and D ha\e been defined as numerical coefficients whose value has to 
be calculated for each medium and for the appropriate wavelength but is always of the order 
of 1. To emphasize the dimensional structure of K and D, we rewrite (65) and (67) in the 
following manner: 

The power-law singularity (64a), entered into the integral (63) yields the build-up factor 
proportional to x^, which was anticipated in (2). This result can be easily derived by replacing 
(Hs—'p)^ with — ^ in the integral (63):* 

B(x)cc~.j'_°'^Jp exp [i,x-p)x]{t.-pt''-' = x'' ^ Jj^ dt exp (-Q(-0"''' - r(g+ 1) ^*^^^ 

(See appendix C, tables 1, 2 for values of K, ixg/i^s, C/jjl,,.) 

If desirable, one can also carry out the build-up factor integration by entering into (63) 

not merely the singularity of y as given by (64a), but the expression of I dsy(p,s,\) as given by 

(61a) or (61b). We may write 

B{x)^ ~ r dp exp [(mo -P)^] (mo- P)-' fl + '^'^^^T, (70) 

and proceed, for example, by expanding the binomial 

r Mo(x-Xo) 7 ^ r(-g+.) r mo(x-Xo) T 

and integrating term by term, which yields 

B{x)ocJ: ^^^^^J [-iio(\-\o)xY=iF,i-K,l-iio[\-\o]x) = L_K{-iio[^-\,]x). (72) 

Here iFi indicates, as usual, the confluent hypergeometric function whose power expansion is 
represented by the 2^, and L-k indicates the Laguerre function, to which the iFi reduces 
when its second parameter equals 1. The binomial expansion (71) is actually unnecessary 
since (70) is directly recognizable as an integral representation of iFi or L-k- Equation (72) 
was derived from a schematized treatment of the straight-ahead problem in (14) [2]. 

The singularit3^ (14b) for \s = ^m, entered into the integral (63), yields a build-up factor 
represented by a transcendental function of unusual type, namely 

I rico 

B(X)CC —^ dp exp {(^rn — p)x + [DfxJ{fXm — p)y^^} 

= £>M«(OM«:r)-2«^ f dt exp {{t+i-"'){Dnr.xy'']=D^„{D^,„xr''^G{[Dix,^xY'\ (73) 

Here t={iJim—p) {Dl^m)~^'^x}'^ , the contour A loops around the origin in the positive direction 
and 

G{2)=^^ f dt exp [{t+t-"')z], (74) 



*The expression of the integral as a r-function is found, for example, in reference [20], last formula on p. 245. 
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For large values of z, the function G can be evaluated by steepest descc^nt method in the 
proximity of the saddle point, which lies at t — 2~'^'^. The result is 



(7(z)-(27r)-^/^(2-^/-l3E)-'/-' exp (2-^/^32). 
This fornnda, entered into (73), yields 



B{,r)oz{Dn,,:){Dy.,nx) 



-5/6<^-5,/6 



(3x)-/2 exp [2-'i'Z{DtXrnfy"] 



(75) 



(76) 



This build-up factor is proportional to x ^^^ exp [H{fjL,nxy^''^], as anticipated in (3), with 

H=2-''''W'''' = ?>{iv'0'l2^,nP^y'\ (77) 

This result was derived in (11 13c) [3] for the straight-ahead approximation, that is with C in 
the place of C . (See appendix C, table 3 for values of ^m, H, D, and C). 

The case \=\s=\m with the singularity (64c), may be treated by taking the second deriv- 
ative of B{;x) with repsect to D^''^, since this operation generates (64c) from (64a), and then 
I'eplacing D wath the appropriate value (68) of /; namely 7^/4, and H with ///2'/^ = 0.63 //. 

The function G{z) has been tabulated throughout the range ^>0 by carrying out tlie 
intc^gration (74) numerically along the steepest descent path lin{t + t~^'^^) = 0. Howevei- no 
application of this tabulation has been required because Spencer [6] found it convenient to 
evaluate intc^grals of the type 



2^^ j dl){^,,-j)Y exp{{fjL,n — p)x + [Dfji,J{iJ^, 



_^.N11/2 1 



-P)V 



(78) 



by expansion into powders of [Dij,„J[/jL,n — p)Y^^ followed by term-by-term integration in the 
manner of (70) and (72). The method proved practical because* the con\ (M"genc(* of the expan- 
sion w^as moderately rapid despite the large vahu* of the argument. (The nund^er of terms in 
the expansion corresponds to the multiplicity of scattering [6]). 

I thank L. V. Spc^nccM- for many helpful discussions and for contributing the appendix C. 
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11. Appendix A. Connection Between the 
X-Ray Distributions From a Plane Iso- 
tropic and From a Point Isotropic 
Source 

The X-ray distribution from a plane isotropic source is 
represented by the solution of eq (5), with /(u^) = 1. The 
equation for a point isotropic monoenergetic source is a partic- 
ular case of (4) in which one takes the source term Sir,u,\) = 
5(r)5(X — Xo). The solution of this equation will be indicated 
as F{r,Ur,\), where Ur=u-rlr is the radial component of u. 

Following up the discussion of plane-symmetrical geom- 
etries in section 1 , we notice that the plane isotropic source 
5(x)5(X— Xo) may be regarded as a plane integral over a 
distribution of point isotropic sources located at ro= (0,?/o,2;o) 



5(a:)5(X-Xo): 



u: 



t(2/„dzo5(r — ru)5(X-\,). (Al) 



US 



Owing to the linearity of the transport eq (4), the distribution 
Y{x,Ux,^) will be similarly expressed as an integral over the 
distribution F{r,Ur,\). Consider an isotropic point source at 
iO,yo,zo) and a point P at (:c,0,0) . The distance from the source 
to P is r = [a:2-f ?/o2+2'o2]i/2 and Ur is equal to (uxX — Uyyo — UzZo)/ 
[x^ + jf/o^ + ^o2]^/2, ' We have 



YixjiA) =j^^j^JyodzoF([x' + yo'+zo'F', 



Uj^c — Uyyo — UzZp \ 

(A2) 



If we take cylindrical coordinates {p,\p,x), with p^ = ?/o^ + ^o^ 
and )/' = arcos [~ iuyyo-\-U2Zo)/p{l — 'u^)^^^, (A2) becomes 

vf.., ^^ r^"^/ f" ^ ;^A^^-L 211/2 UxX+jl-ulY^'pcosxP \ 



(A3) 



The distribution F{r,Ur,\) may be represented by the 
Legendre polynomial expansion 

F(r,Ur,\) = i:ii2m)i4^)-^Fi{r,\)Pi{Ur), (A4) 

and similarly 

Y(x,Ux,\) = j:i(2l+l)i^7r)-^Yi{x,X)Pi(ux). (A5) 

Since the f "" d xpFiHu^x^ (l-n^)^ p cos iA]/[a:2 + p2]^) equals 
simply 2TrPiiux)Pi{xl[x'-^p'']), (A3) reduces to 

Yi(x,X) = jj pdpFi{[x^-\-p^'^,\)Pi{x/[x'^ + p^>^) 

= {" RdRFi{R,\)Pi{xlR). ' (A6) 

The Legendre coefficients with Z = represent the total 
intensity distribution at each point, irrespective of direction 
of propagation. Since Po=l, the relationship (A6) simplifies 
for Z = to 

Fo {x, A) = f " RdR Fo {R, X) ( A7) 

Fo{x,\) =x-^dYo{x,\)/dx. (A8) 

12. Appendix B. Note on the Treatment of 
Strongly Nonself -ad joint Systems 



Consider the system of equations. 

n 
^ G'nmXm \ pXn+ 1 ^^ C^ . 



(Bl) 



In the limit p = this system has a trivially simple solution 
since the first equation contains the variable xi alone and all 
other equations can be solved chainwise. If p is sufficiently 
small, the terms pxn+i may be regarded as small perturbations 
with respect to the trival case p = 0. 

Our method of proceeding is this: If one assumes a trial 
value of the ratio X2lxi, the first equation determines the 
absolute value of xi and X2. Thereafter the n-th equation 
determines the n+l-th variable, chainwise. Owing to the 
smallness of p, any error in the initial estimate of X2/X1 which 

causes a wrong estimate of 2 cinmXm, causes a much larger 

m= 1 

error in the estimate of Xn+i, and so on. A trial and error 
method with successive estimates of X2/x^ may thus proceed 
quite rapidly, because the amplification of the error in 
successive equations soon causes the trial solution to become 
grossly unrealistic, as judged by some suitable criterion. 

Experimental calculations were made with the homoge- 
neous system 



Xn=(Xn-i-\-Xn-2 + Xn-i-Xn+l)/2, 



(B2) 



assuming Xn = for negative n. This system may be regarded 
as representing a game of chance of the parchesee type, in 
which case x must be positive. Successive variables are 
linear functions of the trial value of :ci/xo with alternate and 
rapidly increasing slope. Calculation up to n = 8 suffices to 
show that the solution becomes unacceptable unless X\/xo 
'-0.325. 

As a physical model for the considerations presented here, 
one may think of a chain of pendulums with asymmetrical, 
or ''coaster", couplings, such that the motion of each pen- 
dulum influences the next pendulum on its right strongly, 
but the one on its left only weakly. If the right-to-left cou- 
pling were nil, then the motion of the whole set would be fully 
determined proceeding from the left to the right. Because 
the right-to-left coupling is weak, one can still determine the 
whole motion easily by approximation methods proceeding 
in the same way from left to right. Any error in the initial 
estimate of the motion of one pendulum implies a large error 
in the corresponding estimate for the pendulums further to 
the right. 

We call this situation ''strongly nonself-adjoint" because it 
is a typical feature of the self-adjoint problem with sym- 
metric coupling that one must consider the whole set of pen- 
dulums simultaneously. Similarly, we would characterize 
the trivial case of a completely one-way coupling as "extreme 
nonself-ad joint" . 

In general terms one might state the situation as follows. 
Self-adjoint systems are represented by symmetrical mat- 
rices. Extremely nonself-ad joint systems are represented by 
matrices all of whose elements vanish on one side of the 
diagonal; these systems have simple solutions. Systems with 
matrices having small-valued elements on one side of the 
diagonal should generally be amenable to easy approximate 
treatment, especially if these elements are confined to a strip 
near the diagonal. 

13. Appendix C* Determination of Wick's 
Eigenvalue by Semiasymptotic Itera- 
tion Procedures 

This method for determining eigenvalues can perhaps best 
be illustrated by proceeding step by step through the solution 
of a specific example. Thus, suppose we want to determine 
the eigenvalue of eq (44) for the case of Pb, that is, for 
C/iJirn = 2.90. To begin with, we backtrack a little and define 
a variable u={C/iJLm)^^^v. In terms of this variable, eq (44) 
becomes 



(CI) 



Next, we multiply this equation by powers of (u'^/2), integrate 
over the u space, and thus write down formally moment 
equations of which the first two are 



{C/pim)Koi-Ki=(C/nm)Ko 
(C/fjLm) K, + K2= (C/m.) {K, + Ko) 



(C2) 



where Kn=il/27r)jdu(u^l2)'^XQM- If we write ^ = C/fjLm and 
insert numbers, these two equations become 






(C3) 



Now the right sides of these equations are the moments of 
a smooth, positive, monotonically decreasing (as ?/— >oo) 
function, which we shall call R{u). We shall assume that a 
reasonable representation for this function is R{u)=be-'^^''^^^K 
From (CI) we then hsiwe xo(u)=U + uV2]-^he-P^^'^/^K 
If we take moments of this quantity, we readily obtain the 
relation 



KilKo = g=^ 



l/(^ 



/(^^) 



(C4) 



*By L. V. Sppncer. 
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where f{x)=xe^\ — Ei{--x)]. We now establish an iteration 
j)r()(*eciure usin^' eq (C3) and ((-4), namely 



9 






£(«+i) = 2.90 ^("^'> 
1 



o{n + l)__ 



gin + l)-^l 



(C5) 



where n refers to n'th approximation. The middle equation 
in (C5) is simply a rewrite of the first of eq (C3), while the 
last equation hi (C5) is obtained by fitting moments of the 
assumed R(u)' to the right sides of eq (C3). 

The function f{x) can be easily calculated from standard 
tables. For interpolation between tabulated values we use 
the relation fix ^5) ^f(x)[l-{-8( 1 -f l/a;)]-5. 

To begin the calculation we make a guess that g{^) = \. 
This yields ^(«) = 1.90, /^(")=0.5, and /30^« = 0.95. Continuhig, 
/(0.95) =0.5863. Thus 



^(1) = 1.3407 
^(1) = 1.5593 

,8(i) = 0.42722 



^(2) = 1.4556 
^(2) = 1.4444 
/3(2) = 0.40723 



^«> = 1.4864 
^(3) = 1.4136 
,8(3) = 0.40219 



^(4) = 1.4931 j7('^)= 1.4951 
^(4) = 1.4069 ^('^)= 1.405 
i3(4) = 0.40111 



/(^(i)^(i)) = 0.5172 /(iS(2)^(2)) =0.49284 /(/3«)^(3)) ^ 0.48632 /(/3(4)^(4)) =0.48481 



The difi'erence between successive ^'s is converging by a 
factor ^3-4 at each step. Consequently, ^(°°) should be, 
within a small fraction of a percent, equal to 1.404. This 
yields C=0.484C. This compares with Wick's variational 
value of C=0.499C. (See eq (51).) 

This type of calculation can be extended easily by taking 
into account more of the moment equations (C2). A calculation 
making use of four moment equations yields an approximate 
eigenvalue in agreement with that of the variational method 
to better than 1 percent. 

The iteration procedure for the case of a monotonically 
increasing attenuation coefficient can be set up in the same 
way. Our original intention was to calculate the eigenvalues 



of tables 2 and 3 by this method using four moment equations 
and starting from a variational method value. It turned out, 
however, that the semiasymptotio calculations agreed to 
wdthin about a percent or better with the variational calcu- 
lations in all parts of the tables 2 and 3 and thus we decided 
to present the numerical values of tlui variational calculation 
results obtained from (;ci (51) and (52). 

The agreement between variational and semiasymptotic 
calculations argues strongly for an accuracy of 1 percent or 
so in the eigenvalue determinations. Howc^ver, in practical 
siuations, K and //depend also upon first and second deriva- 
tives, respectively, of fx. These are difficult to determine 
accurately. There may thus be quite sizeal:)l(i errors in the 
tabulated values of //and /) (tal)l(i 3) and of iis/f^s (table 1). 



Table 1. Values of ixslt^s cind C//jls for various materials and 
energies 



Source 
energy 


H2O 


Al 


Fe 


m«/m» 


C//X. 


A»/m» 


C/fls 


A«/m« 


Clus 


Mev 

10 
8 
6 
4 
3 
2 

1 
.8 

'.I 

.3 
.2 
.15 


7.5 
6.6 
5.5 
4.0 
3.1 
2.1 

.96 
.73 
,52 
.32 

.22 
.13 
.086 


7.54 
6.87 
5.96 
4.91 
4.22 
3.38 

2.37 
2.12 

1.86 
1.58 

1.41 
1.23 
1.12 


4.1 
4.0 
3.8 
3.3 
2.7 
2.0 

. 96 
.73 
.52 
.32 

.22 

.155 

.145 


6.21 
5.91 
5.43 
4.68 
4.10 
3.33 

2.37 
2.12 

1.86 
1.57 

1.40 
1.21 
1.09 


"6.'48" 
1.2 
2.0 
2.1 
1.8 

.89 
.72 
.54 
.37 

.32 
.31 
.33 


"4.' 64" 
4.56 
4.25 
3.89 
3.30 

2.36 
2.11 
1.84 
1.52 

1.32 
1.01 
0.767 


Source 
energy 


Sn 


Pb 


u 


10 
8 
6 
4 
3 

2 

1 

.8 
.6 
.4 

.3 
.2 
.15 


"0^40" 
.86 

1.4 
1.1 

.92 
.83 

.77 

.74 
. 72 
!70 


"3." 61" 
3.49 

3.13 
2.22 
1.95 
1.63 
1.17 

.798 
.417 
.223 


0.33" 

1.7 
1.7 
1.6 
1.5 
1.35 

1.2 

.95 

. 77 


2." 88" 

2.65 
1.75 
1.43 
1. 033 
.579 

.333 
.133 
. 0647 


"0.'22" 

2.0 
1.9 
1.8 
1.6 
1.4 

1.2 
.97 

.82 


2.' 72" 

2.48 

1.53 

1.21 
.852 
.449 

.260 
.100 
.0485 
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Table 2. Values of the exponent K=CliJ,s of the build-up factor for various combinations 

of the parameters jlslfJ^s Cind Cf/Jis 

[Values of As/m« and C/fXs for various materials and energies are given in table 1] 



m«/m« 


C/m. = 0.02 


0.05 


0.1 


0.25 


0.5 


1 


2 


4 


6 

189.5 
76.5 
38.5 
15.69 
8.16 

4.24 

2.28 

1.173 

.870 

.671 


8 

272 

108.9 
54.7 
22.2 
11.32 

5.87 
3.11 
1.674 
1.165 
.897 


0.02 
.05 
.1 
.25 
.5 

1 
2 

4 
6 

8 


0. 0200 
. 01939 
. 01846 
.01617 
. 01341 

. 01003 


0. 0516 
.0494 
.0473 
.0411 
.0338 

,0252 


0. 1095 
.1052 
.0988 
.0844 
.0687 

.0506 
.0335 


0.374 
.332 
.293 
.231 
.1798 

.1289 
.0841 


2.63 
1.328 

.871 
.547 
.388 

.265 
.1696 


13.00 
5.49 
2.98 
1.439 

.887 

.559 
.345 


42.3 
17.21 
8.84 
3.81 
2.11 

1.211 
.709 


113.0 

45.5 

23.0 
9.48 
4.96 

2.67 
1.477 
.821 
.577 
.447 











Table S.— Values of fx^ H=S[w'0/2fim(d^fildX^)nY^\ 
D = 4H^/27 and C for various materials 





















C 










CcmVg per 


Material 


fim 


H 


D 


Compton 

wavelength 

unit) 




cmyg 








H20 


0. 0167 


2.0 


L3 


0.118 


AI 


.0216 


2.1 


1.4 


.0942 


Fe 


.0300 


2.8 


3.4 


.0826 


Sn 


.0351 


2.6 


2.7 


.0689 


W 


.0391 


2.5 


2.4 


.0619 


Pb 


.0410 


2.3 


1.8 


.0594 


U 


.0425 


2.1 


1.4 


.0565 
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